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The effect of a Coulombic coupling on the dynamics of a quantum dot hybridized to leads is 
determined. The calculation treats the interaction between charge fluctuations on the dot and the 
dynamically generated image charge in the leads. A formally exact solution is presented for a dot 
coupled to a Luttinger liquid and an approximate solution, equivalent to treating the lead dynamics 
within a random phase approximation, is given for a dot coupled to a two- or three-dimensional 
metallic lead. The leading divergences arising from the long-ranged Coulomb interaction are found 
to cancel, so that in the two- and three-dimensional cases the quantum-dot dynamics is equivalent to 
that obtained by neglecting both the dot-lead Coulomb coupling and the Coulomb renormalization 
of the lead electrons, while in the one-dimensional case the dot-lead mixing is enhanced relative to 
the noninteracting case. Explicit results are given for the short-time dynamics. 

PACS numbers; 71.10.Pm, 73.21.-2b, 73.20.Mf 



I. INTRODUCTION 

The quantum dot, a system comprised of a small 
number of spatially localized levels coupled to one or 
more metallic leads, is one of the paradigmatic prob- 
lems of condensed matter theory and is highly relevant to 
nanoscience. Theoretical studies of quantum dots typi- 
cally involve two competing effects: local interactions, 
which constrain the possible electronic configurations of 
the quantum dot, and hybridization with the leads, which 
mixes different dot eigenstates. In this paper we study 
a third crucial and physically relevant interaction: the 
Coulomb coupling between the charge on the dot and the 
charge on the leads. This interaction is most important 
in the case of a dot weakly coupled to leads, because in 
this case the physical situation is of infrequent transitions 
between states of well defined integer charge. 

Short-ranged dot-lead interactions play a crucial role 
in X-ray absorption lineshapesi and the Kondo effect^ 
and were studied using bosonization methods by Schotte 
and SchottCfS But the general issue of dot-lead inter- 
actions received relatively little attention in the recent 
nanoscience literature, although Gefen and co-workers 
have noted that introducing an additional Coloumbically 
coupled lead into a standard quantum-dot problem may 
lead to a non-Fermi-liquid stated."— 

The interaction we wish to treat may be written as 



Hconi = nd^ j d^rpa{r) 



e 



(1) 



Here we denote by rid the operator giving the number 
of electrons on the dot and by Pa(f) the operator giv- 
ing the number density of electrons at position r relative 
to the impurity in channel a of a Z?-dimensional lead, e 
is a background dielectric constant. Consistency requires 
that a theory involving HqouI involves also considerations 
of the Coulomb efi^ects on the density-density response 
of the lead electrons. In the case of a one-dimensional 



lead, the standard techniques of bosonization^ allow us 
to include both effects, expressing the low-energy elec- 
tronic physics entirely in terms of density and spin fluc- 
tuation operators and enabling a complete and formally 
exact theory. We find that the interaction has (as ex- 
pected from previous results obtained for the X-ray edge 
problemii^) a profound effect on the dynamics of the im- 
purity, qualitatively changing the power laws describing 
the time evolution. We also present an extension to the 
case of higher-dimensional leads. We employ a canoni- 
cal transformation first introduced by Bohm and PineS)^ 
which can only be carried out approximately but which 
captures the essential physics of screening. In this case 
we find that the Coulombic renormalizations of the dot- 
lead coupling and the lead-electron dynamics cancel, so 
the dot-lead problem can be treated within an essentially 
noninteracting electron approximation. 

The physics of local defects in Coulomb Luttinger 
liquids has received previous attention. Fabrizio et 
alJ^ and Maurey and Giamarchiii have considered the 
effects of one or more impurities (described as a short- 
ranged potential scatterer) on a Luttinger liquid with 
Coulomb interaction. Lii>i^ has investigated the screen- 
ing of a test charge, however, in a model with only short- 
ranged electron-electron interactions. None of these au- 
thors treated an impurity with dynamical charge fluctu- 
ations or a Coulombic dot-lead coupling. We also note 
that Lerner et alJ^ have studied a model related to the 
specific model we consider, namely, an impurity adja- 
cent to a Luttinger liquid. They did not consider the 
Coulombic coupling but did include charge fluctuations 
and found a nontrivial structure of the transmission co- 
efficient. 

The effect of a dot-lead interaction on the population 
of a quantum dot has been studied by Goldstein et alt^ 
using density matrix renormalization group and classical 
Monte Carlo simulations. Particular attention was paid 
to the screened case where the long-range part of the 
Coulomb interaction can be neglected. The electronic 
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tunneling was investigated using a Coulomb-gas analysis, 
which allows for an expansion to all orders in the dot- 
lead hybridization. The enhancement of the electronic 
tunneling at the resonance found in Ref. [13 has a similar 
physical interpretation as the results presented in this 
paper. 

The paper is organized as follows. In Sec.|ll]we intro- 
duce the model. Section Hill presents the analysis of the 
one-dimensional problem, including the explicit forms ap- 
propriate to a dot with Luttinger liquid leads. Section llVl 
presents the extension of the results to the case of two- 
and three-dimensional leads. An application of the re- 
sults to the short-time impurity dynamics is presented in 
Sec. IVl We summarize our results in Sec. IVTl 



II. MODEL 

We consider a quantum dot that is hybridized to one 
or more leads, but neglect dot-lead potential scattering 
and do not consider the case where the dot breaks the 
lead into two semi-infinite leads connected only by hy- 
bridization through the dot. These simplifications allow 
us to focus on the consequences of the Coulombic cou- 
pling HqouI, Eq. dl]). The effects we have neglected in- 
troduce additional complicationsiSrJ^ii^ whose interplay 
with the Coulombic coupling will be treated in a future 
paper. 

The quantum-dot problem is described by a Hamilto- 
nian of the general form 



H = H\ 



lead 



dot 



Coul 



(2) 



with -ffcoui given by Eq. ([T|). The quantum-dot Hamilto- 
nian i/dot may be written as 



^^dot = £d nd + ^nd{nd - 1) + 



(3) 



Here U is the dot charging energy and Ud = J2aa '^aa-d'aa 
is the operator giving the total number of electrons on the 
dot. The operator dj^^ creates an electron with energy 
Ed and spin a in dot state a. The ellipsis denotes other 
on-dot interactions, for example, the Hund's coupling J. 

We label the lead orbitals by a and write the lead 
Hamiltonian as 



aha 

1 r „2 



2 J €\r~r'\ 



The dot-lead hybridization is given by 



akaa 



(5) 



The standard quantum-dot physics arises because 
[-ffmix, ^^dot] 7^ so that thc interaction-induced con- 
straints on the dot occupancy interact nontrivially with 
the hybridization, giving rise, for example, to the Kondo 
effect. The physics we wish to investigate arises because 
[Hnux, -ffcoui] 7^ so that a hybridization event changes 
the local charge, giving rise to a long-ranged electric field 
to which the lead electrons react. 



III. ONE DIMENSION 

If the leads are one-dimensional and the electron dis- 
persion is linear, then the low-energy physics of iJicad 
may be expressed in terms of boson o^'^^ and the results 
used to solve the model exactly. We illustrate the method 
here for a single lead with multiple channels; the gener- 
alization to multiple leads is straightforward but involves 
more complicated algebra. 

We imagine a system with linear dimension L (which 
we will later take to infinity) and periodic boundary con- 
ditions so the allowed values of q are 2T:n/L with n ^ 
an integer. We combine spin and orbital quantum num- 
bers into a superindex /3 = 1, . . . , M. The physics is con- 
veniently represented in terms of right (A = and left 
(A = — ) moving particle-hole pairs p^(q), which obey the 
commutation relation [p^ (q), {—q')] = ±Sqq>qL/2Tri^ 
These can be recombined into boson operators. 



n^(9) 



Ppii^ ^ Pp 



which obey the volume commutation relation 

TT 

The total particle density in lead /3 is given by 

Pi3{q) = iq<l>i3{q)- 



(6) 
(7) 

(8) 
(9) 



The lead electron creation operator ipxp may also be 
expressed in terms of bosons as 



(10) 



where c^^^ creates an electron with momentum k and 
spin cr in state a, p{r) is the operator giving the charge 
density at position r, e is a background dielectric con- 
stant and the ellipsis denotes any additional short-ranged 
interactions. It is important that the same long-ranged 
interaction appears in Eqs. ^ and (jH). 



Here wc have omitted the Klein factors that carry the 
Fermi statistics and have introduced a small positive in- 
finitesimal factor rj arising from the correct normal or- 
dering of the operators 1^ 

A key result of the theory of one-dimensional conduc- 
tors is that in thc absence of Umklapp scattering the 
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low-energy physics of the leads may be described by new 
boson operators (f)h, Ilf, related by a linear transformation 
to the operators 0^,11^ and also obeying the canonical 
commutation relations Eq. (|5]). In terms of the new op- 
erators the lead Hamiltonian becomes 



&=1,...,M q 



Kt{q)Uh{-q)nt{q) 



KM 



(t>b{-q)(t)b{q) 



(11) 



with velocity parameters Vb and interaction parameters 
Kb determined by the bare velocities and interactions of 
the lead eigenstates. 

In the most general case the transformation relating 
the 4>b^^b to the 0^,11^ is complicated; in particular the 
equations for (pb and 11^ each may involve both (jjp and 
Tip and the final h combinations need not have a simple 
interpretation in terms of the original densities i/'^'i/'/S- 
However in the most relevant case, where all of the chan- 
nels in a given lead have the same bare velocity and the 
interactions conserve the total lead density, then one of 
the channels (which we take to be & = 1 for definiteness) 
is (up to an overall factor) the total charge density and 
is given by 



4>b=i{q) ^ -^X0/3(g), 
Ilb=i{q) = -l=^n^(q). 



Conversely, for any of the original indices j3 we have 

1 



(l^piq) 
mq) 



M 
1 



06=1 (g) + 

Tib=l{q) + 



(12) 
(13) 

(14) 
(15) 



with the ellipses representing the other operators (all 
commuting with 4>b=i, IIb=i) needed to make up the full 
operator. In particular, the electron annihilation opera- 
tor assumes the form^ 



V'A/3(a;) = 
with 



g '/ml 



/27r?7 



X e 



(17) 



with C and D the transformation coefficients which di- 
agonalize the Luttinger-liquid Hamiltonian. 

Comparison to Eq. ^ shows that the total charge den- 
sity p is related to (f)b=i by 



(18) 



p(q) = VMiq(l3i{q). 



Thus writing 



In 1 



q' 



(19) 



(20) 



with A the inverse of a short-distance cutoff we find that 
the long-ranged Coulomb interaction between conduction 
electrons in the Luttinger liquid is 



MnvpVc 2 



:q^Mq)M-q)Wg. (21) 



Here we introduced a dimensionless measure of the 
Coulomb interaction strength 

V, = . (22) 

In a general Coulomb-coupled Luttinger liquid we have 
also a short-ranged part of the interaction, parametrized 
by dimensionless constants 171^2 such that if the Coulomb 
interaction were negligible we would have 



Ki,sitiq) = 



1 + .91 - 92 



(23) 



1 + .91 + ff2 ' 

Vi,sR{q) = VFKi,sR{q){l + 91+ 92) ■ (24) 
Including the Coulomb interaction gives 
Ki,sB{q) 



Ki{q) 



1 



1+91 +32 



vi[q) = Vi^sR{q)\ 1 + — ■ . 

V 1 + .91 + 52 



(25) 



(26) 



The dot-lead interaction, Eq. ([T]), is transcribed into 
the new representation as 



H, 



TTVp 



Coul 



AIL 



Y,iqMq)MV,Wgnd. (27) 



(28) 



The linear coupling between the dot occupancy Ud and 
the lead density in Eq. (f27l) may be removed by a 
canonical transformation which shifts <j)i{q) — (t>i{q) — 

= ^-lI^MV.W, ^ . 

' vi{q) " 1 + gi + 92 + MVcWg 

In the shifted Hamiltonian the dot-lead interaction is 
eliminated, i/ioad retains the form of Eq. ([TT|) and the 
dot energy and the local interaction U are decreased 
by A and 2A respectively, with the polaron shift A given 
by 

2M L ^ VFKi(q) « 

q 



2ML Y 1 + 51 + 52 + MV.Wq 



(29) 
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Here the polaron shift gives the static interaction be- 
tween the dot charge and the image charge it induces in 
the lead. 

The canonical transformation acts on an operator O 
by C ^ e^-^Ce"'-^. From Eq. © we sec that 



S 



-nd 



^Vni(-<z)^. 
ML ^ n 



(30) 



Under the canonical transformation the fermion operator 
ip\j3{x), Eq. ([TC)). becomes 



while the operator dj,^. becomes 



V'A/3(a;) 



4a, e" 



(31) 



(32) 



The factor multiplying ip is purely imaginary. In this 
paper we will only need to consider combinations d^ipip^d 
for which this factor and a similar one coming from the 
commutator needed to combine the factor multiplying 
d with the boson operators in cancel. We may thus 
write the transformed operator appearing in the dot-lead 
hybridization term as 



[•*• 



ni(9) 



(33) 



Thus in a dot Coulombically coupled to a one- 
dimensional conductor the explicit Coulombic dot-lead 
coupling may be eliminated by a canonical transforma- 
tion. The physics associated with the Coulomb interac- 
tion is expressed via a renormalization of the dot-lead 
hybridization. 

As will be seen below, for the evaluation of physical 
quantities the crucial objects are lead-operator expecta- 
tion values of the form 



F(x,t) = (6/3(2;, 0^^(0,0)), 



(34) 



where is the right-hand side of Eq. p3p with the d^ 
operator removed. 

For later use we present explicit formulae for the two 
most interesting cases: a metallic nanotube lead and a 
single-channel Luttinger liquid with SU(2) spin rotation 
invariance, specializing further to the case of a local dot- 
lead hybridization (x = 0). In both of these cases the in- 
teractions and velocities for the noncharge channels may 
be approximated by the free-fermion values of vt = vp 
and Kh = 1. It is convenient to multiply and divide 
by the free-fermion correlator Fq, which is proportional 
to 1/t at long times and has an appropriate short-time 
cutoff. We obtain 



F(i)=Fo(i)e*(*\ 



(35) 
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u - 



V 
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FIG. 1: (Color online) Real and imaginary part of the renor- 
malization factor <I>(t), Eq. (|36p . for a Luttinger liquid with 
M = 4 channels and Coulomb cutoff (tube diameter) A — 
10 /rj plotted against time (in units of bare cutoff rj/vp). Solid 
blue line: Re$(t) computed from exact expression with di- 
mensionless coupling 14 — 0.9. Dashed green line: Re3'(t) 
computed from Eq. (|36p but with screening factor Zq set to 
zero. Dotted red line: Im${t) computed from exact expres- 
sion with Vc — 0.9. Dash-dotted purple line: Im<I>(t) com- 
puted with screening factor Zq set to zero. 



where 



2n ^ 

TTT. 2^ 



1 - e-'""!* -Bq{l- e 



9#o 



2|g| 



and 



q = «l('7)kl, 



Bq = 



(1 



VF\q\, 

ry \2 



(36) 

(37) 
(38) 



In these formulae the Luttinger-liquid correlations are 
expressed by the factors Ki{q) and ujg/uj';^ while the effect 
of the dot-lead interaction is carried by Zg. 

The real and imaginary parts of $ are plotted in Fig.[T] 
for parameters appropriate to a nanotube (see below). 
In the absence of the Coulombic dot-lead interaction 
{Zq = 0), one finds Re $(t) < reflecting the suppression 
of the electronic tunneling by interaction effects in one di- 
mension. However, Eq. shows that the Coulombic 
dot-lead interaction acts to reduce the magnitude of the 
negative term in $. We see from Fig. [T] that the dot-lead 
coupling in fact changes the sign of Rc$(t), therefore 
enhancing the dot-lead hybridization above the nonin- 
teracting value. The long-time asymptotic behavior of 
Re(j){t) is \og{vpt/r])/M in the screened case {Zq ^ 0). 

A metallic nanotube has two conducting channels and 
two spin states, so M = 4. For a metallic nanotube 
with e = 1 and bare Fermi velocity vp = 5.3 eV A, 
one has Vc ~ 0.9. We also remark that the parame- 
ter A in Eq. ([^0]) is on the order of the inverse of the 
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tube diameter, which is much greater than the basic lat- 
tice constant. Further, in metalhc nanotubes the short- 
ranged interactions (and indeed all interaction effects ex- 
cept the long-ranged Coulomb interaction) are believed 
to be negligiblei^ so that Kc^sR = 1 and Vc,SR. = vp and 



nanotubc 



(<?) 



1 4- MVcWn 



(39) 



.nanotubc 



(g) = VF^Jl + MVWq- 



Finally, evaluating the polaron shift we find (choosing 
the cutoff A to be the inverse of the nanotube diameter 

^nanotubc) 



An 



tubo[cV] 



5.6 

'^nanotubc [-^] 



(41) 



For a typical nanotubc diameter of 12 A, one obtains 



■^nanotubc 



0.5 eV. 



IV. LEAD DIMENSION GREATER THAN 1 

The treatment of the previous section relied on a par- 
ticular feature of one-dimensional (non-nested) systems, 
namely that the electronic degrees of freedom could be 
entirely eliminated in favor of a set of effectively noninter- 
acting boson excitations, one of which is the density. In 
dimension higher than one a complete elimination of elec- 
tronic degrees of freedom is not possible but a separation 
of density fluctuation and fermionic variables may be ef- 
fected via a canonical transformation method introduced 
by Bohm and Pines.— The canonical transformation can- 
not be carried out exactly but an approximate imple- 
mentation can be performed, which has the same level 
of accuracy as the familiar random phase approximation 
(RPA). This approximate implementation leads to a the- 
ory which is very similar in form to the one-dimensional 
theory derived above, but of course expressing the char- 
acteristic physics of higher-dimensional systems. 

The analysis proceeds from the Hamiltonian in Eq. © . 
Following Bohm and Pines^ we introduce canonically 
conjugate boson fields Pq and Qq obeying (in D spatial 
dimensions) 



(42) 



To begin we consider states \'ip) which are a direct product 
of fermion and boson eigenstates and restrict attention 
to states satisfying 



PqlV) = 0. 



(43) 



We then may shift the density operators of all lead states 
a, 



Pa(q)^pf'*(q)^Pa(q) + 



(44) 



where N denotes the number of lead states and de- 
notes the Fourier transform of the Coulomb interaction 
47reV(g^i^) in £) = 3 and 2e^TT / {\q\L'^) in £> = 2 dimen- 
sions. Provided that we consider only wave functions \tp) 
which obey the subsidiary condition Eq. (|43)) the Hamil- 
tonian written in terms of p^^^ift jg equivalent to the orig- 
inal Hamiltonian. 

Bohm and Pines now introduce a canonical transfor- 
mation to shift Pq by ^^/V^X^a /°a(q)- This transforma- 



(40) tion is effected hy H ^ e^^He'''^ with 



(45) 



aq 



After this transformation, the subsidiary condition for 
the wavefunctions becomes [Pq — ^/^X^a Pa(q)]l^') ^ 0: 
which allows us to replace X^a /'^(q) by Pq/-\/V^. 

Under the transformation, i/dot remains invariant 
while iJniix becomes 



E 



(46) 



and the excitonic dot-lead coupling HqouI is 

i^coui = J2 \^Pci^d- (47) 



Similarly, the lead-electron Green function Gioad(r,t) = 
(c,,(r,t)cL(0,0)) with Caa{r,t) = J d^ge"i--c„q<,(t) be- 
comes 



Glead(r,0 Gioad(l',i) 



(48) 



with Gioad computed with the transformed Hamiltonian. 
Bohm and Pines show that the renormalization implied 
by Eq. (^5)) is in essence the RPA reduction in the elec- 
tronic spectral weight. 

These are exact results. As far as is known, the trans- 
formation of the remainder of the lead Hamiltonian can 
only be carried out approximately, for example, by ex- 
panding the exponentials to obtain a series of multiple 
commutators. Keeping the exact first-order commutator 
and approximating the second-order term by its vacuum 
expectation value, hq, Bohm and Pines obtain 



i/load ^ fkclkaCakCT - /^Qqq ' Jq (49) 

akcT aq 

q 

with the electron current operator jq given by 

*q • jq = E^'^k+q - ek)clkaCa(k+q)a (50) 



kcr 
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and the plasma frequency flp defined in terms of the elec- 
tron stress-energy tensor by 



X! '^o (*^k+q + Ck-q - 2ek) 



(51) 



The coupling between Qq and the divergence of the 
fermion current, and the terms dropped in the approxi- 
mate canonical transformation are to be treated pertur- 
batively. Neglecting them is equivalent to treating the 
plasmon as an undamped boson and retaining the term 
in the electron self energy which comes from the electron- 
plasmon coupling. Adding the iQq(q-jq) term in leading 
order of perturbation theory restores the plasmon damp- 
ing found within the RPA while the additional neglected 
terms give beyond-RPA physics. 

These transformations have reduced the problem to 
one analogous to that solved in the previous section. 
We now shift the field Pq by —^/Vqnd to remove the 
dot-lead coupling. The canonical transformation O — > 
e^^Oe^"^^ with S ~ j-u X]q v'^Qq'^d precisely cancels 
the Coulomb- induced renormalization of the lead fermion 
operator, so that in the transformed variables the dot- 
lead hybridization takes the unrenormalized form 



(52) 



Thus in dimensions greater than 1 and within the RPA 
approximation the Coulombic dot-lead coupling compen- 
sates for the Coulomb-induced reduction in electronic 
spectral weight, so that Coulombic effects drop out of 
the tunneling problem (except for the 'image charge' or 
polaron-shift reduction in the dot energy and screening 
of the dot interaction). The tunneling effects may be cal- 
culated using free electrons (that is to say, computing the 
local lead correlators without including the RPA or GW 
self energy) . 



V. QUANTUM-DOT DYNAMICS 

The previous sections have shown how to reduce the 
Coulombically coupled dot-lead problem to an expansion 
in powers of the dot-lead hybridization which is appropri- 
ately renormalized by bosons. The lowest-order term in 
this expansion implies a master-equation approach to the 
dot dynamics, which we use to illustrate our formalism. 

For simplicity, we consider a quantum dot with a sin- 
gle (non-degenerate) level described by a density matrix 
which is diagonal in the occupancy basis and is written 
as 



Pd = Po 



+ A|1)(1|. 



(53) 



Here Pq and Pi denote the occupation probabilities of the 
empty state |0) and the singly-charged state Insert- 
ing Eq. (|53p into the von Neumann equation of motion 
for the density matrix, expanding to leading nontrivial 



order in the hybridization and using Pi = 1 — Pq yields 
the master equation 

Po(t) = / dt' ( Ri^a{t ~ t') 



Po{t') [Po^i(t - t') + i?i^o(< - <')] (54) 



with the transition probabilities 



i?o->i(t) = 2|rpRe [P(t), 
i?i_o(t) = 2|rpRe [F{t)e''''] , 



(55) 
(56) 



where F{t) is defined in Eq. (f34|) . 

We note in passing that because here we restrict at- 
tention to thermal equilibrium, the long-time limit of the 
populations is a steady state determined by the detailed- 
balance condition, i.e. Pq (as a function of £d/T) is a 
Fermi distribution in the limit t — > oo. Of course, at zero 
temperature T = 0, orthogonality effects may cause the 
system not to equilibrate. 

We first examine these equations in the Markovian 
limit in which the dot dynamics are slow enough and the 
kernels F decays fast enough that Po(t) may be treated 
as a constant and extracted from the integral in Eq. (|54p . 
This yields 

Po{t) = 7^l^o - Po{t) (7^o^l + 7^l^o) , (57) 



where 

7^o^l(erf) =2|rpRc 



dTP(r) ( 



(58) 



and TZi-fo{ed) — TZo-fii—ed) are the corresponding 
Golden-Rule transition rates. 

To obtain an idea of the effects of screening we approxi- 
mate the logarithmic functions tUq and K{q) by constants 
oj and K. The integrals may then be performed analyti- 
cally and we find at zero temperature 



for times t ^ rj/v with 



M 



(59) 

(60) 
(61) 



if the Coulombic dot-lead interaction is included, and 



Y 



2M 



K-\-K~ 

a = K 2M 



(62) 
(63) 



if not. 
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FIG. 2: (Color online) Time evolution of the occupation prob- 
ability Po{t) for the symmetric case, Ed = 0, with initial values 
Po(0) = 0.51 and Pi(0) = 0.49. Solid line (blue online): Po{t) 
computed from Eq. (|54p with dimensionless coupling 14 — 0.9. 
Dashed line (green online): Po{t) computed from Eq. (|54[) but 



with screening factor Zg set to zero. We assume a tunneling 
amplitude T = 0.1 vp/rj (corresponding to a bare tunneling 
time To = 1677/?'^) and the same parameters as chosen in 

Fig.m 



Because for repulsive interactions < if < 1 we see 
that in the presence of Coulombic dot-lead coupling we 
have Y < I, whereas if the coupling is neglected we have 

y > 1. 

The rates are then determined as 



l£rf| 
v/t] 



Y-1 



a 0{-£d) 
T{Y) ■ 



(64) 



Equation ((M]) is of the form of a basic tunneling rate, 



1 

To 



(65) 



times a factor expressing the effect of correlations. We 
see that if the Coulombic dot-lead coupling is neglected, 
the interactions suppress the tunneling rate, whereas in 
the presence of Coulombic dot-lead coupling the relax- 
ation rate is enhanced. In the symmetric case Ed ~ 
the Markov rate vanishes if the coupling is neglected but 
diverges if is retained. The steady state of the system 
obtained from the master equation, Eq. (j57p . depends on 
the value of Y. In the most interesting case, Ed = 0, 
it is known that although the Markov rate vanishes for 
y > 1, for y < 2 it remains the case that as t —J- 00 
the occupancy Po(t) tends to the thermal-equilibrium 
value, Po{t) 1/2, while for y > 2 the system does 
not equilibrate at zero temperaturei^ In the nanotube 
case of interest here the effective exponent K{q) vanishes 
as g — > so that at sufficiently long scales the model 
without Coulombic dot-lead coupling would fail to equi- 
librate; however because the increase is only logarithmic, 
for reasonable nanotube parameters the effective expo- 
nent would only become greater than 2 for q ^ A. 



Figure [2] presents the time evolution of -Po(t) obtained 
by solving the master equation, Eq. ([M)) . for a dot 
coupled to a nanotube with and without the dot-lead 
Coulombic coupling. We have considered the particle- 
hole symmetric case Ed = for which the equilibrium 
value is P{Q) = 1/2 and have begun the simulation in 
a non-equilibrium initial condition. The strong enhance- 
ment of relaxation by the Coulombic dot-lead coupling is 
evident. 



VI. CONCLUSIONS 

In summary, we have studied the effect of a Coulom- 
bic dot-lead coupling on the dynamics of a quantum dot. 
This coupling is always present, but its effects seem not 
heretofore to have been examined. We find that it has 
an important effect on the dot-lead dynamics. Two cases 
have been considered: a dot coupled to a Luttinger liquid 
and a dot coupled to two- or three-dimensional metal- 
lic leads. The effects are particularly profound for one- 
dimensional leads. It is well known from previous work 
that in the absence of the Coulombic dot-lead coupling 
the Luttinger liquid correlations of a one-dimensional 
lead strongly reduce the dot-lead hybridization, leading 
(at low enough scales and for strong enough interactions) 
to a complete suppression of tunneling and failure of the 
system to equilibrate at T = 0. The dot-lead Coulomb 
interaction is shown to overcompensate for this effect, 
leading to a divergence in the dot-lead hybridization. 

The enhancement of the electronic tunneling due to the 
Coulombic dot-lead interaction has the following quali- 
tative interpretation. The suppression in the usual case 
comes because when adding a charge one has to push 
the other charges in the Luttinger liquid aside, while the 
presence of the screening interaction implies that other 
charges need not be pushed away because a screening 
cloud has to be formed. 

We presented estimates for parameters appropriate to 
a carbon nanotube, which is one of the most widely used 
one-dimensional leads. We also presented explicit for- 
mulae which could be used for more detailed numeri- 
cal simulations along the lines of Refs. [13 and [U For 
three-dimensional leads, the effects are found to be less 
dramatic, but still significant: the consequence of the 
Coulombic dot-lead coupling is that the Coulombic renor- 
malizations drop out of the problem, so that the dot con- 
ductance should be studied using lead Green functions 
unrenormalized by the RPA or GW corrections to the 
electron effective mass and scattering. 

Our results rely on several approximations. The most 
crucial is that density fiuctuations in the leads can be 
represented as noninteracting bosons. The standard re- 
sults of Luttinger liquid theory^ justify this approxima- 
tion for the case of one-dimensional leads (at least in the 
universal low-energy limit) while for the case of higher- 
dimensional leads our approximations are at the same 
level as the random phase approximation. 
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Our paper leaves several avenues for future research. 
Our explicit results are perturbative in the dot-lead hy- 
bridization. A numerical or analytical treatment to all 
orders, leading in particular to an expression for the lin- 
ear response I — V curve, would be very valuable. An 
extension of the work to the non-equilibrium case of non- 
zero bias voltage is also important. 
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